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We present the microscopic theory of improper multiferroicity in BiMnOa , which can be summa- 
rized as follows: (1) the ferroelectric polarization is driven by the hidden antiferromagnetic order in 
the otherwise centrosymmetric C2/c structure; (2) the relativistic spin-orbit interaction is respon- 
sible for the canted spin ferromagnetism. Our analysis is supported by numerical calculations of 
electronic polarization using Berry's phase formalism, which was applied to the low-energy model 
of BiMnOs derived from the first-principles calculations. We explicitly show how the electric po- 
larization can be controlled by the magnetic field and argue that BiMnOa is a rare and potentially 
interesting material where ferroelectricity can indeed coexist and interplay with the ferromagnetism. 



Introduction. Today, the term 'multiferroics' is typ- 
ically understood in a broad sense, as the systems ex- 
hibiting spontaneous electric polarization and any type 
of magnetic ordering.^ Such materials have a great po- 
tential for practical applications in magnetic memories, 
logic, and magnetoelectric sensors, and therefore at- 
tracted enormous attention recently. Beside practical 
motivations, there is a strong fundamental interest in 
unveiling the microscopic mechanism of coupling be- 
tween electric polarization and magnetic degrees of free- 
dom. Nevertheless, the combination of ferroelectricity 
and /erromagnetism, what the term 'multiferroicity' was 
originally introduced for, is rare. Such a combination 
would, for example, provide an easy way for manipulat- 
ing the electric polarization P by the external magnetic 
field, which is directly coupled to the net ferromagnetic 
(FM) moment, etc. The canonical example, where spon- 
taneous electric polarization was believed to coexist with 
the FM ground state, is BiMnOs. However, the origin 
of such coexistence is largely unknown. Originally, the 
ferroelectric (FE) behavior in BiMnOa was attributed to 
the highly distorted perovskite structure stabilized by the 
Bi6s "lone pairs''.^ However, more resent experimental 
studies (Ref. |3|) and first-principles calculations (Ref. U) 
suggested that the atomic displacements alone result in 
the centrosymmetric C2/c structure, which is incompat- 
ible with the ferroelectricity. In our previous papers 
(Refs. la and [g) we put forward the idea that the ferro- 
electricity in BiMnOa could be improper and associated 
with some hidden antiferromagnetic (AFM) order. The 
purpose of this work is to provide the complete quantita- 
tive explanation for the appearance and behavior of the 
FE polarization in BiMnOa. 

Method. The basic idea of our approach is to construct 
an effective Hubbard-type model 
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for the MnSd-bands near the Fermi level and to include 
the effect of all other ("inactive") states to the defi- 
nition of the model parameters of the Hamiltonian %. 
Thus, the model is constructed in the basis of 40 Wan- 



nier functions in each unit cell (including three t2g- and 
two Cg-orbitals for each spin and for each of the four 
Mn-sites), by starting from the electronic structure in 
the local-density approximation (LDA). The Greek sym- 
bols denote the combination of spin and orbital indices. 
All parameters of % are defined rigorously, on the basis 
of the density functional theory (DFT). The details can 
be found in the review article (Ref. \J\j and in our pre- 
vious papers (Refs. [^ andQ). Briefly, the one-electron 
part {&) is derived by using the generalized downfold- 
ing method. One of important parameters in i"^ is the 
large (about 1.5 eV) crystal-field splitting between two 
eg-levels, which is caused by the Jahn- Teller distortion 
and manifests itself in the orbital ordering. The screened 
Coulomb interactions iUaji-ys) are obtained by combining 
the constrained DFT technique with the random-phase 
approximation (RPA):^ namely, the screening by outer 
electrons (such as 4sp-electrons of transition metals) and 
the change of spacial extension of the atomic wavefunc- 
tions upon the change of occupation numbers can be eas- 
ily taken into account by solving Kohn-Sham equations 
within constrained DFT approach. On the other hand, 
the "self-screening" by the same type of electrons, which 
contribute to other bands due to the hybridization ef- 
fects (for example, the Sd-electrons in the oxygen band 
will strongly screen the Coulomb interactions in the 3d- 
band near the Fermi level) , is included in the perturbative 
RPA treatment. The self-screening is very important in 
solids and substantially reduces the value of the effective 
Coulomb repulsion U (defined as the screened Slater in- 
tegral F°) in the 3d-band of manganites.^ In BiMnOa, it 
is only about 2.3 eV,— that has important consequences 
on the behavior of interatomic magnetic interactions. 

The model ([1]) is solved in the Hartree-Fock (HE) 
approximationii 



^k + V) |C„k) — fnklC'nk), 



where £k is the Fourier image of £ij = ||i^ || and, if neces- 
sary, includes the relativistic spin-orbit coupling (SOC), 
V is the self-consistent HE potential, and |C„k) is the 
eigenvector in the basis of Wannier functions (where the 



plane A' 'v 

1 



plane .v'z' 



^ ^^ ^ 

I2 A 1 4 

f o mf^ 

4 Aw!2 S^ 



3/-)^ 3r--)- 




Ir-r' 33'--r^ ^rfsr^ 



FIG. 1. (Color online) Schematic view on the orbital order- 
ing and corresponding interatomic magnetic interactions in 
the pseudocubic x'y' and z'x' planes. In the unit cell of 
BiMnOs, there are four Mn sites (indicated by numbers), 
which form two inequivalent groups: (1,2) and (3,4). The 
nearest-neighbor FM interactions Jjviv operate in the hatched 
bonds. The atoms involved in the long-range AFM interac- 
tions Jlr are denoted by arrows. The inversion centers are 
marked by '*'. 



spin indices are included in the definition of n)^ 

Once the orbital degeneracy is lifted by the strong lat- 
tice distortion, the HF theory provides a good approx- 
imation for the ground-state properties. The effect of 
correlation interactions, vi^hich can be treated as a per- 
turbation to the HF solution,- on the magnetic ground 
state of manganites is partially compensated by the mag- 
netic polarization of the oxygen states: if the former tend 
to stabilize AFM structures, the latter favors the FM 
alignment. - Due to this compensation, the mean-field HF 
theory, formulated for the minimal 3d-model, appears to 
be rather successful for the ground state of manganites. 
Magnetism and the inversion symmetry breaking. 
First, let us explain the main idea of our work.^'^ What 
is the possible origin of multiferroic behavior of BiMnOa 
and how can it be controlled by the magnetic field? 

(1) The lattice distortion leads the orbital ordering, 
which is schematically shown in Fig. [1] in two pseudocu- 
bic planes (the orbital ordering in the y'z'-plane is similar 
to the one in the z'x'-plane). This orbital ordering pre- 
determines the behavior of interatomic magnetic interac- 
tions, which obey some general principles, applicable for 
manganites with both monoclinic (C2/c) and orthorhom- 
bic (Pnma) structure;^!^ namely: besides conventional 
nearest-neighbor interactions (shown by hatched lines), 
one can expect some longer-range (LR) interactions be- 
tween remote Mn-atoms, which operate via intermediate 
Mn-sites. These sites are shown by arrows. 

(2) Why should the LR-interactions exist? The answer 
is directly related to the fact that the on-site Coulomb 
repulsion U is not particularly large. Therefore, besides 
conventional superexchange (SE), there are other inter- 
actions, which formally appear in the higher orders of 
the 1/ [/-expansion and connect more remote sites. This 
mechanism is rather similar to the SE interaction via in- 
termediate oxygen sites, but the role of the oxygen states 



is played by the unoccupied eg-orbitals of the intermedi- 
ate Mn-sites.^ By mapping HF total energies onto the 
Heisenberg model, one can obtain the following param- 
eters of interatomic magnetic interactions:^'^- Jnn ^ 5 
and 6 meV (where slightly different values correspond to 
different bonds) and Jlr '^ —3 nieV. Thus, these inter- 
actions are at least comparable. Besides them, there are 
finite (of the order —1 meV) interactions in the bonds 1-2 
and 4-4 across the inversion center, which finally define 
the type of the magnetic ground state of BiMnOa. 

(3) Without SOC, the LR interactions tend to stabilize 
the AFM tut structure (where the arrows denote the di- 
rections of spins for the four Mn-sites in the unit cell). 
This AFM order destroys the inversion centers (shown 
by '*' in Fig. [T]) and thus should give rise to the FE 
polarization. Since the tilt structure satisfies the sym- 
metry operation T<Si{my\'R3/2} (where my is the mirror 
reflection y -^ —y associated with the one half of the 
monoclinic translation R3, and T in the nonrelativistic 
case flips the directions of spins, which are not affected 
by ruy), P is expected to lie in the zx-plane.^^ 

(4) Thus, the FE behavior in BiMnOa should be caused 
by the AFM order. However, this conclusion seems to 
contradict to the FM ground state of BiMnOs.^, The 
contradiction can be reconciled by considering the rel- 
ativistic SOC, which is responsible for the weak ferro- 
magnetism. Since the FM component is additionally 
stabilized by the isotropic interactions Jnn, the ferro- 
magnetism is not so "weak" , and the resulting magnetic 
structure, obtained in the HF calculations for the low- 
energy model, is strongly noncoUinear (Fig. [5]). It be- 
longs to the space group Cc, where the only nontrivial 
symmetry operation is {toj,|R3/2} and the magnetic mo- 
ments in the relativistic case are transformed by rriy as 
auxiliary vectors. Thus, the net FM moment is aligned 
along the y-axis, while the x- and z-components form 
the AFM structure. Other magnetic configurations have 
higher energies. The details can be found in Ref. 0. 

By summarizing this part, the C2/c symmetry of 
BiMnOa is spontaneously broken by the hidden AFM 
order. The true magnetic ground-state of BiMnOa is 
strongly noncoUinear, where the FM order along the y- 
axis coexists with the AFM order, and related to it FE 
polarization, along the x- and z-axes. Our scenario not 
only explains the coexistence of ferroelectricity and ferro- 
magnetism, but also shows how the electric polarization 
P (and the symmetry of BiMnOa) can be controlled by 
the external magnetic field 'B—{Q,By,0) coupled to the 
FM moment. This basic idea was formulated in Ref. % 
In the present work we are able to provide the numerical 
estimates for P and to discuss its behavior in details. 

Electric polarization. Since the crystal structure of 
BiMnOa has the inversion symmetry, there will be no 
ionic contribution to P, and the main mechanism, which 
will be considered below, is of purely electronic origin. 
In principle, the magneto-elastic interactions in the tiit 
structure may cause the atomic displacements away from 
the centrosymmetric positions and give rise to the ionic 
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FIG. 2. (Color online) Fragment of crystal and magnetic 
structure corresponding to the lowest HF energy. The Bi- 
atoms and indicated by the big light grey (yellow) spheres, 
the Mn-atoms are indicated by the medium grey (red) spheres, 
and the oxygen atoms are indicated by the small grey (green) 
spheres. The directions of spin magnetic moments are shown 
by arrows. The inversion center is marked by the symbol '*'. 
The left lower part of the figure explains the orientation of 
the Cartesian coordinate frame. 



term. Nevertheless, such calculations would require the 
full structure optimization, which cannot be easily in- 
corporated in the model analysis. The first-principles 
calculations for HoMnOa show that electronic and ionic 
terms are at least comparable. ^^ Thus, we expect that the 
electronic contribution alone could provide a good semi- 
quantitative estimate for P. Moreover, the behavior of 
electronic contribution presents a fundamental interest as 
it allows one to explain how P in improper multiferroics 
is induced solely by the magnetic symmetry breaking. 

The modern theory of electric polarization allows one 
to relate the change of P to the Berry phase of Bloch 
eigenstates»i^"— It is particularly convenient to use the 
formulation by Resta, where the Berry phase is com- 
puted on the discrete grid of k-points, generated by the 
N1XN2X N3 divisions of the reciprocal lattice vectors 
{Ga}.J^ Then, the position of each point in the Brillouin 
zone is specified by three integer indices {0 < Sa < Na). 

k.......3 = ^Gi + — G2 + — G3, 

and components of the electric polarization in the curvi- 
linear coordinate frame formed by Gi, G2 and G3 can 



be obtained as^ 
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and similar expressions hold for 72 and 73. Eq. ^ im- 
plies that the only meaningful quantity in the bulk is the 
polarization difference between two states that can be 
connected by an adiabatic switching processJ^"— 

In the present case, S = ||(C„k|C'n'k')ll is the overlap 
matrix, constructed from the HF eigenvectors |C„k) in 
the occupied part of spectra, taken in two neighboring 
k-points: k^ks^^s^^sj and k'^ks^+i^s^^^g for 71, etcJ^ 
The polarization ([2]) was first computed in the curvilinear 
coordinate frame and then transformed to the cartesian 
frame shown in Fig. [2P^ In all the calculations, we used 
the mesh of 72 x 72 x 36 points in the Brillouin zone. 

Without SOC, the AFM alignment of spins at the sites 
1 and 2 yields finite polarization. However, the symme- 
try of the system also depends on the magnetic config- 
uration in the sublattice 3-4. As discussed above, the 
electric polarization in the tiit structure lies in the zx- 
plane {P^ = 2.1 ^C/cm^ and Pz = 0.1 /iC/cm^). The 
tiit structure can be transformed to the titi one with 
the same energy by the symmetry operation {C^|R3/2} 
(where Cy is the 180° rotation around the y-axis), which 
changes the direction of P: P^ 
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other hand, the t4-4-4- structure (which has higher en- 
ergy) is transformed to itself by {C^|R3/2}, and cor- 
responding electric polarization will be parallel to the 
2/-axis {Py = 4.8 /iC/cm-^). Other magnetic structures, 
characterized by the FM alignment of spins at the sites 
1 and 2 (such as tttt £^nd ttti)! preserve the inversion 
symmetry and result in zero net polarization. 

Without SOC, one can easily evaluate separate contri- 
butions to P of the states with the spin t and |. For 
the tiit structure, the vector of the electric polarization 
takes the following form: P''^-'- = ^{Pj;,±Py,Pz), where 
Py = 5.7 fiC/cm^, and the values of Px and Py are listed 
above. This result is very natural, because the distribu- 
tion of the electron density for each spin does not have 
any symmetry and, therefore, the electric polarization 
P'''"'' has all three components. On the other hand, the 
electron density with the spin f in the tJ-it AFM state 
can be transformed to the one with the spin J, by the sym- 
metry operation {TOy|R3/2} and, therefore, Pj = —Py- 
Thus, in the total polarization P = P^ -I- P^, the x- and 
z-components with different spins will sum up, while the 
largest y-components will cancel each other. 

Furthermore, one can evaluate the individual contri- 
butions to P coming from the i2g-band, which is sepa- 
rated by an energy gap from the eg-band.— This yields: 
P^^o = -0.8 ^C/cm2 and P^^" = -0.3 ^iC/cm^. Thus, 
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FIG. 3. (Color online) Magnetic-field dependence of tfie elec- 
tric polarization, the angle </> between spin magnetic moments 
at the Mn-sites 1 and 2, and the vector of magnetic moment 
at the site 1 (shown in the inset). 



the i2g-band is polarized opposite to the Cg-band, that 
substantially reduces the value of P. 

The SOC results in the canting of spins away from the 
coUinear tiit state and towards the FM configuration. 
It will reduce the value of P. In the HF ground-state (see 
Fig. [2]), the angle cj) between spin magnetic moments at 
the sites 1 and 2 is 137°, and the electric polarization is 
reduced till P^. = 1.6 /xC/cm^ and P^ = 0.1 ^C/cm^. 
This effect can be further controlled by the magnetic 
field, which is applied along the y-axis and saturates the 
FM magnetization. Since the absolute value of the local 
magnetic moment is nearly conserved, the increase of the 
FM component along the y-axis will be compensated by 
the decrease of two AFM components along the x- and 



z-axes. The corresponding FE polarization will also de- 
crease. Results of HF calculations in the magnetic field 
are shown in Fig. [31^^ Sufficiently large magnetic field 
(^35 Tesla) will align the magnetic moments at the sites 
1 and 2 ferromagnetically (0=0) and restore the C2/c 
symmetryjS The electric polarization follows the change 
of and complete disappears when 0=0. However, the 
decline of P is much steeper: for example, Px and Pz 
are reduced by factor two already in the moderate field 
By ^ 5 Tesla, corresponding to ^ 100°. Moreover, P^ 
is always substantially smaller than P^. 

Concluding remarks. We have proposed the micro- 
scopic theory of improper multiferroicity in BiMnOs, 
which is based on the inversion symmetry breaking by 
the hidden AFM order. We have estimated the FE po- 
larization and explicitly shown how it can be controlled 
by the magnetic field. Our scenario still needs to be 
checked experimentally, and apparently one important 
question here is how to separate the intrinsic ferroelec- 
tricity in BiMnOa from extrinsic effects, caused by the 
defects. For example, the values of the FE polarization 
obtained in the present work, although comparable with 
those calculated for other improper ferroelectrics on the 
basis of manganites,^^ are substantially larger than the 
experimental value 0.062 ^C/cm^ (at 87 K), which was 
reported so far for BiMnOa.— Nevertheless, we believe 
that systematic study of manganites with the monoclinic 
C2/c symmetry and finding conditions, which would lead 
to the practical realization of scenario proposed in our 
work, presents a very important direction, because it 
gives a possibility for combining and intermanipulating 
the /erroelectricity and /erromagnetism in one sample. 
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